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Abstract 

We develop the pQCD description of diffraction excitation of heavy flavours in DIS and we 
derive the analytic formulas for the mass spectrum in leading logm^. The result illustrates 
nicely non-factorization properties of the QCD pomeron. We predict a very steep rise of 
the charm content of diffraction dissociation of photons at small xjp. We evaluate the con- 
tribution of open charm to scaling violations in the structure function of the pomeron. 
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Following Ingelman and Schlein [IJ , Regge factorization |2j is often applied to diffraction 
dissociation (DD) of (virtual) photons 7* + p — > X + p' into states X of mass M (large 
rapidity gap (LRG) events), which are reinterpreted as a deep inelastic scattering (DIS) on 
pomerons radiated by the target proton, endowing the pomeron with the usual attributes 
of a particle such as the partonic structure function F 2 tp(P,Q 2 ) and the flux of pomerons 
<j>Tp{xjp) / Xjp in the proton (fll], [|: 



(M2 + g2 ^(7*-X) 



dt dM 2 



■bw{x w )F 2 w{p,Q ) ■ (1) 



t=Q 167T Q 2 

Here Q 2 is the virtuality of the photon, W and M are c.m.s. energy in the photon-proton 
and photon-pomeron collision, (3 = Q 2 /{Q 2 + M 2 ) is the Bjorken variable for the lepton- 
pomeron DIS and xjp = (Q 2 + M 2 )/(Q 2 + W 2 ) = x/ (3 is interpreted as the fraction of the 
momentum of the proton carried away by the pomeron. 

The Ingelman-Schlein model has never been derived from a QCD analysis. Quite to the 
contrary, the unequivocal conclusion from the QCD approach to DD is the non-factorization 
of the QCD pomeron jJJ: the xjp dependence of dap can not be reabsorbed entirely in the 
Q 2 , f3 and flavour independent pomeron flux function 0ip(xip), whereas the (3, Q 2 and flavour 
dependence of dan can not be contained entirely in the // structure function of the pomeron. 
Different aspects of the non-factorization in DD have been discussed in [|, §, |, 0, ||; 
the non-factorizable colour dipole approach to DD 0, 0| is well known to provide very 
good quantitative description of the HERA data on LRG events f|, [Toll . Regarding the 



applicability of pQCD, the crucial observation J| is that in DD of transverse photons, the 
qq pairs have a, Q 2 -fndependent, typical transverse size ~ 1/rrif. This allows to quantify 
the factorization breaking in DD into heavy flavours on a more quantitative basis than for 
the (predominantly nonperturbative) DD into light flavours. The demonstration of this 
factorization breaking is the subject of the present paper. 

One of the main points of the present communication is a derivation of the f3 dependence 
of the pQCD factorization scale, 

This f3, rrif dependence of the pQCD factorization scale drives the flavour and (3 dependence 
of the flux of pQCD pomerons in the proton and nicely illustrates the breaking of Regge 



factorization for the QCD pomeron: each and every flavour brings in a new, and explicitly 
/3-dependent, flux function into the menagerie of pomeron fluxes in the proton. We predict a 
very steep rise of the charm abundance in DD towards small xip, which is a distinctive feature 
of the QCD approach to DD as compared to the Regge factorization models. We present 
also an evaluation of the contribution of open charm excitation to the counterintuitive rise 
with Q 2 of the pomeron structure function at large values of (3. Finally, we comment on the 
even more dramatic factorization breaking in DD of longitudinal photons, where the pQCD 
factorization scale q$ ~ \Q 2 entails the Q 2 -dependent "flux of pomerons" which can be 
related to the Q 2 dependence of the gluon structure function of the proton (for the similar 
situation in DD into dijets see [El). 

DD into open heavy flavour, X = qq, is described by diagrams of Fig. 1. It dominates 
at M 2 ~ Q 2 considered here. The relevant formalism is described in detail in j|, |6|. The 
mass of the state X is given by M 2 = (m) + k 2 )/z(l - z) , where m/ is the quark mass, k 
is the transverse momentum of the quark with respect to the 7*-pomeron collision axis and 
z is the fraction of light-cone momentum of the photon carried by the (anti) quark. Other 



useful kinematical variables are e 2 = z(l — z)Q 2 



m 2 f and 



k 2 +^ = (k 2 +m}) 



M 2 + Q 2 
M 2 



(3) 



After the standard leading logK 2 resummation, the cross sections of the forward (t = 0) 
DD of (T) transverse and (L) longitudinal photons takes the compact form 0. M 



dM 2 dk 2 dt 



t=o 



—e f a em a s {q ) 



m 2 + k 2 



M 3 cos 9 J M 2 -4m 2 



k 2 + m 2 



f 



$ 2 + m 2 <!> 2 2 , (4) 



da j 



dM 2 dk 2 dt 



t=Q 



—e f a em Q a s {q J 
o 



(mj + k 2 ) 3 



M 7 cos 9 J M 2 -4m 2 



$ 2 . 



(5) 



Here 9 is the quark production angle with respect to the 7*-pomeron collision axis, 



k 



k 



+ 



2kK 2 



k 2 + e 2 Va 2 - b 2 a 2 - b 2 + aVa 2 - b 2 



$2 = [^fixjp,^ 2 



Va 2 - b 2 k 2 + e 2 



(6) 
(7) 
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a = e 2 + k 2 + k, 2 , b = 2kn and f(xjp, k 2 ) = dG(xjp, k 2 ) / d\ogn 2 is the unintegrated gluon 
structure function of the target proton. In the derivation of the pQCD factorization scale 
we follow the analysis ||: At small k 2 < q 2 the expression in the square brackets in the 
integrand of (|6|) equals 2ke 2 K 2 /q 6 and tends to a constant value at k 2 > q 2 . Then, in (Q) 
and ([FD one has a logarithmic k 2 integration with q 2 being the upper limit of integration. 
Consequently, q 2 emerges as the pQCD factorization scale (it has already been used as such 
in the running strong coupling as(q 2 ) in (§) and (|)) and to the leading logg 2 , 



2kM A [Q 2 {k 2 + m 2 ) + M 



2 mV 



$1 = (Q 2 + M^(k 2 + m 2 r G(XlP ' q ] (8) 



^ ^ -r „ t/A ,» y J *,, bf M j 2 

$ 2- 77^ I iJ^WTU I _2vj G(s!P,gj. (9) 



M 4 [(£; 2 + m 2 )(M 2 - Q 2 ) - 2m}M 2 
(Q 2 + M 2 f(k 2 + m 2 

Notice a zero of the da L at (k 2 + m 2 f )(M 2 - Q 2 ) = 2m 2 f M 2 . Equations (|), (|) and 
exhaust the derivation of the mass spectrum. 

Subtleties of DD in QCD are clearly seen already for real photons (Q 2 = 0), where 



dM 2 dk 2 dt 



1 2 2 2 / 2,7 2\^-y2/ 2 , ; 2\ 2 

= -7T e^a em a s [m^ + k )G (xjp, + k )mj 



x 



om 2 : 8m 4 - 



' (10) 



M 3 {k 2 + m 2 ) 4 cos JM 2 - Am) 



Evidently, the fc 2 -integrated daT/dm 2 is dominated by the contribution from k 2 < m 2 . 
Consequently, to the leading logm 2 one can take the fc 2 -independent factorization scale 
go = iTij- (The dominance of of the contribution from k 2 < m 2 in ot holds also at large Q 2 , 
entailing the pQCD factorization scale q 2 of Eq. (|2|)). We find 



dan 



dM 2 dt 



VA f~* m2f , if M 2 -Am 2 <m 2 
n"e f a em a s {m f )& (x w ,m } ) x \ ■ ( n ) 

I idfwi if M2 » 4 ™/ 

The M 2 and xip dependence in ([11]) do factor and one may try to reinterpret the xjp 
dependent factor as the pomeron flux function 

ASfit \ fG(xTP,m 2 f )\ 2 

W) = (gk^J ' (12) 

subject to the normalization 0jp (20 = 0.03) = 1 @. The explicit flavour dependence in 
(JT^) breaks the Regge factorization, hence the universal flux of QCD pomerons does not 
exist. 



Consider now the DIS regime of large Q 2 . For transverse photons, the expression in the 
curly braces in Eq. (|4]) is a function of only j3 and Q 2 and can be a basis for the definition 
of the // pomeron structure function F-jp (/3, Q 2 ). The factor a|(go) = a s{ m }{^- ~ fi)^ 1 ) 
also is a function of (3 only and as such can be reabsorbed into the F^\(3,Q 2 ). The xip 
dependence comes entirely from G 2 (xjp,ql) = G 2 (xjp,m 2 (l — and here by virtue of 

the QCD scaling violations the xjp and (3 dependences are inextricably entangled, which 
breaks the Regge factorization (HD explicitly. One can think of a generalized factorization 
at best, leaving the pomeron flux factor to depend explicitly on f3 and on flavour: 

(//) _ ( G{x*,m){l-p)-i) \* 

^ {x ^ P) -[ G ( Xo , m 2 (l-f3)-i)) ■ (13) 

Because q 2 rises with /3, the pQCD is applicable better at smaller M; for the specific case 
of exclusive vector meson production see |TT| . 

At this point, a brief digression on the still more striking, and different, Regge factor- 
ization breaking in the longitudinal cross section is in order [12|. Substituting (|8|) into (|5|), 
one readily finds 



da L 



dM 2 dk 2 dt 



= -TT 2 eia em Q 2 a 2 s (k 2 + e 2 )G 2 (x, k + e 2 ) 
t=o 6 

M[(k 2 + m 2 )(M 2 - Q 2 ) - 2m)M 2 ] 2 
X cos 9{Q 2 + M 2 f{k 2 + m)f^M 2 - Am) 

which decreases with k 2 only as k~ 2 . While <jt is dominated by the contribution from 
k 2 < m 2 , for ol the dominant contribution comes from large k 2 ~ \M 2 — m 2 . This has 
two major implications: First, <tl has the higher twist Q 2 dependence <tl oc <tt/Q 2 oc 1/Q a 
p. Second, in ol the factorization scale ql ~ \Q 2 1 the xjp and Q 2 dependence in daL are 
inextricably entangled and the generalized flux of pomerons does explicitly depend on Q 2 , 
^ipt^ip) G 2 (xjp, \Q 2 ), rather than on (3 and flavour in the case of dar- 

The exceptional case is the triple pomeron region of (3 1 dominated by DD into the 
qqg.. states, where the conditions of the Regge factorization are fulfilled at Q 2 > 3 GeV 2 
for all the flavours simultaneously and the corresponding flux function fjp(xjp) is flavour 
blind. This conclusion readily follows from an analysis of the colour dipole content of the 
triple-pomeron coupling in ref. [13 . 



Hereafter we focus on DD of transverse photons into cc states, which dominates the 
open charm excitation in LRG events at f3 >0. 1-0.2. For open charm the factorization 
scale (0) is still not large, we evaluate the open charm cross section in the colour dipole 
gBFKL formalism described in detail in ||, [5], [Tj. In Fig. 2 we show our results for the 
total cross section of diffraction excitation of cc pairs. The evaluation of this cross section 
requires the t-integration; one can argue that for the heavy flavour excitation the diffraction 
slope B is smaller than the light flavours one (B e i ~ lOGeV -2 ). In the following we use 
B c5 = 6GeV~ 2 (for instance, see The closer analysis of Eqs. (|) suggests a simple 

interpolation ot oc l/(Q 2 + 4m 2 ) between the real photoproduction (Q 2 = 0) and DIS 
(Q 2 3> 4m 2 ) at fixed value of the natural variable x Pom = (Q 2 + 4m 2 ) /{W 2 + Q 2 ). This 
scaling is demonstrated in Fig. 2 where we plot (Q 2 + 4m 2 )ojT c ' ) as a function of the Bjorken 
variable x and xp om . 

Notice the very steep rise, ot oc Xp € om , at xp om < 1CT 3 , with the exponent e ~ 0.72 
which is very close to the asymptotic gBFKL prediction e = 2Ajp = 0.8 (for the prediction 
of the precocious onset of the BFKL behaviour in the charm structure function of the proton 
see ||15|). In Fig. 3 we show the flux function (pjp^ixjp) defined for the above integrated cc 
excitation cross section at large Q 2 , which is dominated by the contribution from (3 ~ 0.5. 
For the region 0.1 < x^p < 10~ 4 accessible at HERA, the convenient parameterization is 

#w )= (^r(fE±i!ir (is) 

\XjpJ \x + p 3 J 

with pi = 0.7233, ]?2 — 0.3939,^3 = 2.377 • 10 -3 , which is different from the flux func- 
tion 0ip(xip) for the valence light-flavour component of the pomeron (pi = 0.569,^2 — 
0.4895, P3 = 0.153 • 10 -3 ) and the flux function /ip(a;ip) for the sea (triple-pomeron) com- 
ponent of the pomeron ( p\ = 0.741, p2 = 0.586, p% = 0.8 • 10 -3 ) 0. The charm content of 
DD is predicted to rise by one order of magnitude from xjp = 0.01 to xjp = 0.0001 (see also 
Fig. 6). 

Now we focus on the Q 2 and (3 dependence of the open charm excitation. The variation 
of the factorization scale (^) at (3 < 0.5 is marginal, but at 1 — (3 <S 1 the factorization- 
breaking (3 dependence of the generalized flux function flTo] ) is quite strong. The distortion 
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factor 

m = (i 6 ) 

<A c W,/3 = 0.5) 

presented in Fig. 4, shows how the shape of the j3 distribution varies with xjp in defiance 
of the Regge factorization. Here we have evaluated T((3) using the GRV gluon structure 
function |I5], for other parameterizations of parton densities the results for the factorization 



breaking will be very similar. 

Despite having discredited the very concept of the pomeron structure function, we can 
not help but use this language to make the contact with what has unfortunately become 
a common presentation of the experimental data on DD. The t- integrated cross section 
measured at HERA can be represented as 

{M +Q ] dM* " 16^ Fd{x ^ & Q } (17) 

with the valence light (qq) and valence charm (cc) decomposition of the non-factorizing 
" diffractive structure function" (we omit the negligible bb contribution, neglect the marginal 
difference between fluxes for the ss and uu, dd excitation and limit ourselves to the transverse 
structure function.) 

F D (x*, (3, Q 2 ) = Q 2 ) + 4&(xr, P)F^(P, Q 2 ) + Mx JP )F^ ea \f3, Q 2 ) . 

(18) 

The result for light flavours, F^ q \(3) w 0.27/3(1 - f3), evaluated with the diffraction 
slope B e i = 10 GeV -2 , provides an excellent description of the experimental data P, pT0f| . In 
Fig. 5 we show our predictions for the (3 dependence of the cc component of F D (xjp, j3, Q 2 ) 
at x w = 0.001. Because M 2 >4m 2 c , the F ( f ) (x JP ,(3,Q 2 ) vanishes at (3 > (3 C = Q 2% m 2 ■ The 
impact of this threshold effect on F^ c \f3,Q 2 ) was for the first time discussed in ||. Fig. 5 
updates Fig. 10 of ref. [Q; in the present calculation we use the more modern gBFKL dipole 



cross section of 17, 11 



In Fig. 6 we show the threshold effect due to opening of the charm production in 
F D (xjp, P, Q 2 ) considered as a function of Q 2 at different values of [3. Here we neglect 
the possible scaling violations in the light flavour contribution Fp(xjp 7 /3, Q 2 ) (see below). 
In the domain of 0.35 < (3 < 0.8 and 5 < Q 2 < 200 GeV 2 of the current HERA experiments, 



the predicted threshold rise of Fd(xjp, (3, Q 2 ) is quite strong and must be observable. The 
error bars of the presently available HERA data || [10[ on Fd(xjp, f3, Q 2 ) are still too large 
for the observation of the charm threshold effect. Fig. 6 clearly demonstrates the rise of 
the charm content of DD towards small xjp: at the typical (3 ~ 0.5 and xjp = 10~ 2 not 
shown here the charm contribution and the threshold effect are ~ 3%, which rises to ~ 9% 
at xjp = 10~ 3 and ~ 25% at xjp = 10~ 4 , cf. different fluxes in Fig. 3. At even smaller xjp, 
not accessible at HERA, the charm content of DD levels off. The strong rcip-dependence of 
the charm content of DD is a non-negotiable consequence of pQCD; such a xip-dependence 
is absent in Regge models [TJ, [L9], p0| . Neither works |[8|, [L9| discuss the strong impact of 
the charm threshold on the Q 2 dependence of DD at large (3. One of the two models of DD 
into open heavy flavour discussed in |2(J assumes a pointlike pomeron-quark coupling; in 
this model too the charm content does not depend on xjp. Such a pointlike pomeron-quark 
coupling is not born out in our pQCD approach. 

Several more comments on our results are in order: 

First, for heavy flavours the results for the (3 distribution are exact in contrast to the 
nonperturbative DD into light flavours, where the j3 dependence is not pQCD calculable. 
Fig. 5 shows that at very large Q 2 3> 4m 2 , in a broad range of /3, the (3 dependence of 
F^) C \xsp 1 (3, Q 2 ) follows the approximation oc (3(1 — (3) fairly well. In terms of the mass 
spectrum, it corresponds to da^f 1 / dM 2 ~ M 2 /(Q 2 + M 2 ) 3 . The departures from this law 
were noticed already in (see Fig. 3 in [[|]). Indeed, in a very narrow region of (3 — > 1, 
i.e., 4m 2 <C M 2 <C Q 2 , from Eqs. (I),©,© one rea dily finds (modulo to the pQCD scaling 
violations) dag/dM 2 ~ M 4 /(Q 2 + M 2 ) 4 and Fg(x W) f3,Q 2 ) oc (1 - (3) 2 . Although this 
observation is hardly of practical significance, the law Fd(xjp, (3, Q 2 ) oc (1 — (3) 2 will be 
applicable also to the light flavour excitation in a very narrow domain for (3 — > 1, in which 
the factorization scale (0) is sufficiently large for the pQCD applicability (for a hint at such 
a (3 dependence see also Fig. 3 in 0]). 

Second, the factorization breaking and the /3-dependent pQCD factorization scale have 
certain implications for the Q 2 -evolution of F^(xjp, (3,Q 2 ), which remains one of the open 
issues in the theory of DD. In || it was shown that an approximate GLDAP evolution 



S 



is recovered in the double leading-logQ 2 , leading- log ^ approximation. Here we wish to 
comment that Eq. (2) implies that DD at large ft is dominated by excitation of qq pairs 
with the transverse size r qq - which decreases with ft, 

ri« -1(1-/9). (19) 

m f 

At large values of the Bjorken variable ft the parton distributions decrease with Q 2 for the 
radiation of gluons. Evidently, the radiation of gluons by a colour singlet system of size r qq - 
requires the condition Q 2 r 2 q 3> 1. Then, Eqs. (PD,(|T9"D hint at the possibility that the larger 
ft the larger Qq oc (1 — ft)^ 1 is needed for the onset of the conventional pattern of the Q 2 
evolution of F§ {x w , ft, Q 2 ) at Q 2 > Q 2 Q . 

Third, the predicted steep rise of the charm excitation cross section has a non-negligible 
impact on the xjp dependence of F D (xjp, ft, Q 2 ). Consider the approximation of this 
quantity in the vicinity of xjp ~ 10~ 3 typical for the current HERA experiments. For 
excitation of light flavours q = u,d,s, our prediction [|7| for the flux function (p-sp(xjp) 
corresponds to e(uds) ~ 0.15, which is close to the soft pomeron model predictions [|[ 
[181, HI, Poll - For the charm excitation we predict e(c) ~ 0.72 at xip < 10 -3 . Although 
the abundance of charm is numerically small, the steeply rising contribution of the charm 
to Fd(xjp, ft, Q 2 ) substantially renormalizes upwards the exponent e for the total DD cross 
section: at ft ~ 0.5 and Q 2 above the charm threshold our crude estimate is e(uds+c) ~ 0.20. 
For comparison, in the triple pomeron region of ft -C 1 our prediction 0] for the flux function 
fjp(xjp) corresponds to e ~ 0.32. Evidently, the estimates for these exponents depend on ft 
and on the range of xjp considered, see Figs. 3 and 4. 

To summarize, we have presented the pQCD derivation of the mass spectrum and cross 
section of DD into open charm in DIS at HERA. Our results unequivocally demonstrate 
strong breaking of the Regge factorization in diffraction dissociation of photons, which 
has already been advocated by two of the authors in a previous work ||. We predict 
a very steep rise of the charm content of DD towards small xjp, which is a distinctive 
consequence of pQCD and can be tested at HERA. We predict strong xip dependence of the 
ft distribution in open charm production at ft — > 1, which also can be tested in the HERA 
experiments. We have found a substantial charm threshold effects in the Q 2 dependence of 
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the diffractive structure function at large (3 which must be observable with higher statistics 
data at HERA. This charm threshold effect leads to a counterintuitive rise of the large-/3 
diffractive structure function with Q 2 , which is stronger at smaller xjp. Our finding of the 
/3-dependent pQCD factorization scale casts shadow on the GLDAP description of the Q 2 - 
evolution of the diffractive structure function at large (3. 
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